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Abstract. We present a closed system of coupled first order differential equations governing 
the secular linear momentum loss of a compact binary due to emitted gravitational waves, with 
the leading order relativistic and spin-orbit perturbations included. In order to close the system, 
the secular evolution equations of the linear momentum derived from the dissipative dynamics 
are supplemented with the secular evolutions of the coupled angular variables, as derived from 
the conservative dynamics. 
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1. Introduction 

The inspiral of compact objects in binary systems is driven by the emitted gravitational radiation 
which carries energy, angular and linear momenta away from the source. Global conservation 
of these quantities implies a radiative orbital evolution, including a possible recoil of the 
system. The energy and magnitude of orbital angular momentum of the binary determine 
the quasi-Keplerian orbit in the post-Newtonian (PN) regime. These quantities are known to 
high accuracy for binaries on eccentric orbits, with PN, spin-orbit (SO), spin-spin (SS), mass 
quadrupole - mass dipole (QM) and magnetic dipole - magnetic dipole (DD) coupling terms pp. 
Moreover, the orientation of the (quasi-precessing) plan of motion, defined by the direction of the 
Newtonian orbital angular momentum, is known to high accuracy with all the above-mentioned 
contributions. 

Due to asymmetries in the configuration of the source the radiation is often emitted 
anisotropically. The linear momentum loss of the binary leads to the recoil of the center of 
mass in the opposite direction, which in extreme situations results the kick-off of the binary 
from its host galaxy. This effect, however, vanishes for equal mass binaries due to the fact that 
the leading order relativistic contribution to the linear momentum loss scales with the mass 
difference dm = m\ — m% of the orbiting bodies. 

The gravitational recoil of a binary system and the classical linear momentum loss of the 
final black hole were discussed in [21 [3] with the inclusion of the lowest multipoles needed for 



the computation of momentum ejection. The first quasi-Keplerian analytic studies were given 
in [3], where the linear momentum flux of gravitational waves from a binary system of two 
point masses in Keplerian orbit were calculated. 1 PN corrections to the gravitational recoil 
were discussed in j5j, while 2PN recoil effects for binaries on quasicircular orbits in [6]. Several 
numerical estimates for the kick velocity of non-spinning binaries were given, e.g. in [U [9], 
with the maximum recoil in this case found as ~ 175 km/s [10] . Recently it has been shown 
[11] , that the ringdown phase acts as an anti-kick as compared to the inspiral and plunge, the 
global result being consistent with the numerical estimates. 

The rotation of the components adds however additional structure to the emitted radiation 
and recoil. Spin contributions to the linear momentum loss are analyzed in [T2], and more 
recently in [T3]. The SO contributions scale with the magnitude of spins which could be high 
for galactic black holes. Numerical analyses indicate that significant gravitational recoil can be 
obtained in spinning binaries [141 \TE\ \16\ \TT\ [TH] even with a high degree of symmetry in the 
configuration, i.e. for equal-mass binaries with antialigned initial spins in the orbital plane [19] . 
In a detailed review [20] . Racine, Buonanno and Kidder have computed the instantaneous linear 
momentum flux emitted by spinning binaries at 2PN order with the inclusion of the next-to- 
leading order SO, SO tail and SS terms. Moreover, the recoil velocity as a function of the orbital 
frequency was given for quasicircular orbits. 

In the present work we give the previously unknown secular expressions for the linear 
momentum loss of an inspiralling binary system, with the inclusion of leading order relativistic 
and SO contributions. The analytic approach presented here is suitable for characterize 
analytically the dependence of the recoil on binary parameters. 

Section 2 contains elements of the conservative dynamics with the inclusion of first post- 
Newtonian and spin-orbit contributions. At the end of the section the secular conservative 
evolution of the relevant angle variables is given. In Section 3 we start from the expressions of 
the SO contributions to the instantaneous linear momentum loss given in [12], then we derive 
the dissipative secular linear momentum evolution. The two sets of secular evolutions couple to 
a closed system of first order differential equations. 

The secular evolutions are derived as follows: first we rewrite the instantaneous evolutions in 
terms of the radial parametrization of quasi-Keplerian orbits [lj. As a result, all the radial 
functional dependencies are expressed in terms of a single variable x, the generalized true 
anomaly. Averaging these expressions over a radial period becomes particularly straightforward 
in terms of the complex version of x, which renders the problem to the computation of residues in 
the origin of the complex parameter plane [21] [22] . This procedure smears out short timescale 
effects and we obtain a simpler dynamics, suitable for monitoring the secular changes. The 
procedure is similar to our previous computations of the SO-induced secular changes of the 
energy, magnitude of orbital angular momentum and relative angles among the orbital angular 
momentum and spins [23] . An additional difficulty in the present computation arises from the 
vectorial character of the linear momentum, and the subsequent system of coupled differential 
equations. 

Our description is generic, being valid for generic (non-circular, non-spherical) orbits. 
Notation. For any vector V its magnitude is denoted as V and its direction (a unit vector) 
as V. 

2. Elements of the conservative dynamics 

In this section we derive useful elements of the conservative dynamics, originating from the SO 
coupling. We derive the result in terms of: 

(a) the physical parameters of the binary: total mass m = mi + 777-2, mass ratio v = m-zjmx < 
1, symmetric mass ratio n = fi/m (where fi = mim^/m- is the reduced mass), dimensionless spin 




(b) dynamical constants of motion: the energy E, the magnitude L of the orbital angular 

1 /2 

momentum L and A = (G 2 m 2 [i 2 + 2EL 2 / 11) (this would be the length of the Laplace-Runge- 
Lenz vector of the Keplerian motion characterized by E and L) 

(c) angular variables related to the orbital momenta: the relative angles of spins among 
themselves 7 and with the orbital angular momentum m, the polar angles ifti of the spins in 
the plane of motion (measured from the intersection 1 of the plane of motion with the plane 
perpendicular to the total angular momentum J = L + S, where S = Si + S2), and 

(d) angular variables characterizing the orbit: the inclination a of the orbital plane with 
respect to the plane perpendicular to J, the angle cp n between 1 and an inertial axis x taken in 
the plane perpendicular to J, finally the angle i/j p span by the periastron An and 1 (see Fig 1). 
These three angles will be referred occasionally as Euler angles, as three consecutive rotations 
with (f> n , a and ip p about the axes z, x and again z transform from an inertial system with J as 
the z -axis to the system with the x-axis pointing towards the periastron and the y-axis in the 
plane of motion. 




Figure 1. The angular variables. 



2.1. Lagrangian dynamics 

The Lagrangian characterizing the dynamics to 1.5 PN orders (tail terms are omitted) is 

C = Cn + £pn + C-so , (1) 



with the contributions 
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Here a = ivSi + v~ 1 S2 is the mass- weighted spin vector. The SO part of the Lagrangian applies 
to the Newton- Wigner spin supplementary condition and was given first in Ref. pQ. The 
Lagrangian gives 

a = a N + a PN + a S o , (3) 



with 
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2.2. Constants of motion 

The constants of motion are the energy E and the total angular momentum J. As the total 
orbital angular momentum undergoes a pure precessional motion |24[ [25] , its magnitude L is 
also a conserved quantity. 

The total orbital angular momentum L = rx (dC/dv) can be decomposed as 



L = Ln + Lpn + Lso 
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with 
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We note here the useful approximate relation 
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2.3. The evolutions of ip and cp n 

The polar and azimuthal angles 9, <p of the reduced mass particle in the inertial system with J as 
the z -axis can be related to the Euler angles (4> n , a, ip) characterizing a non-inertial reference 
system with the x-axis at the location of the reduced mass particle and the y-axis in the plane 
of motion |26j : 

sin 9 cos (p = cos a sin (p n sin ip + cos <p n cos ip , (14) 
sin 9 sin ip = cos a cos (p n sin ip — sin <p n cos ip , (15) 
cos 9 = sin a sin ?/; . (16) 

Tedious but straightforward algebra leads to 

sin a cos ip ■ ■ 

ip = cos a ip — 9 + cos a (p n , 
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cot a ■ tan ib d , 
sm 9 cos ip sm a dt 



The Newtonian orbital angular momentum, expressed in terms of (9, ip) gives 

L 2 N = /j 1 2 r i (9 2 + sin 2 9 ip 2 ) , (L N ) 2 = fir 2 sin 2 9 ip . (18) 
By employing Eqs. (fTT|) . (fT6|) and (Ln) 2 = -Latcoso we obtain the evolutions: 

sin 9 9 = k sin a cos ip ( 1 PN n ~to SO ) , (19) 

sm 6* ip = — j cos a 1 — ^ . (20) 
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The minus sign in the first equation was chosen after taking the square root in such a way that 
when we insert these equations into Eqs. (fT7|) . we obtain 

^ = ~2 f 1 ~ APA 2^ 2 ASQ ) + cos a ' ( 21 ) 
• tanip d 
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sm a at 

which reduce to the correct equations ip = L//j,r 2 and (p n =const at Newtonian order. 

Provided that we can obtain c\ by a complementary method, the evolutions of ip and (p n can 
be given explicitly. This will be done in the following subsection. 



2.4- The evolution of the orbital inclination 
From cos a = J • Ln, as J is conserved, we get 



|r(cosa) = J-^(Ln) • (23) 



dt v ' dt 

Employing Eq. ([3]) gives 
i N = Fxa=^(2-r ) )rL N + ^{(^r-v)[r.(4S + 3a)) - f (4S + 3a) } . (24) 



From here we can derive the time derivative of Ln • We then find that the leading order evolution 
of a is of order S 2 /JL: 
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(We have employed = L, valid to leading order accuracy.) The scalar products can be 
rewritten in terms of the angular variables as 

r • Si = rSi sin «j cos(^ p + X~^i) > 
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In Eq. (|25p . as all terms are of S 2 / J order, it is allowed to employ the Newtonian true anomaly 
parametrization of the orbit: 
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With these expressions, all the desired angular evolutions can be given explicitly as functions of 
X (which, to Newtonian order is ip — if) p , see Figure 1). 

2. 5. The secular evolution of the Euler angles 

We conclude this section by giving the secular evolutions of the Euler angles. The secular 
evolution of any function / (x) is defined as < / >= T -1 J 27r / (x) X 1 d>X- As the instantaneous 
angular evolutions (|22p and (|25p have no Newtonian contributions, we can employ the Newtonian 
period of the orbital motion T = 2irGm[[i/(—2E)] 3 / 2 in calculating their orbital average, 
obtaining: 
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We introduced here the shorthand notations Aip = ip2~ ip\ and tp = (ipi + 1P2) /2. 

We also compute the change in the angle ip over a radial period by the above indicated 
method, as Aip = f? n ip (x) X" 1 By employing Eq. ()2ip and the relevant termd3 form Eq. 



1 The contributions (dt/dx) PN and (dt/dx) so are corrected by a global sign and the factor fl/L, respectively. 



(18) of Ref. PQ we get: 
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with 
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The leading order term gives J Q 27r (L/fir 2 ) (dt/dx) N dx = 2-7T. Defining / ijj p \ = (A^ — 2ir) JT 
we get the secular periastron precession rate as 
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Eqs. ([28]) and ([32]) together with the evolution of the angles ki given by the first two Eqs. (2.17) 
of [23j and the projections of J = L + S1 + S2, determining tpi in terms of the other variables form 
a closed system of differential equations for the conservative evolution of the angular variables. 

3. Elements of the dissipative dynamics 

The instantaneous loss in the linear momentum due to the gravitational radiation can be 

(4) (3) (3) (3) 

expressed as Pj = —{-^lijkljk + ^^ijkljiJki}, where Ijj., 1^ and Jm are the mass quadrupole, 
mass octupole and current quadrupole moments |12j . The secular loss in the linear momentum 
is: 
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with the leading order radiative (N) and leading order radiative SO contributions given by 
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4. Concluding Remarks 

We have derived secular evolution equations, Eqs. (j28|) and (|32p for the angular variables 
a, (fi n and ifi v characterizing the orientation of the orbit, with the inclusion of the leading order 
relativistic and spin-orbit coupling contributions. Together with the evolution of the angles Ki 
given by the first two Eqs. (2.17) of [23J and the projections of J = L + Si + S2, determining 
tfii in terms of the other variables, these form a closed system of differential equations for the 
conservative evolution of the angular variables. 

We have complemented these with the equations expressing the respective secular losses of the 
linear momentum components, Eqs. (|33p . Then we have a closed system of coupled first order 
differential equations, giving the linear momentum loss during the inspiral in analytic form. 

This system of equations is also suitable for numerical evolution, which in principle can lead 
to the value of the recoil velocity right before the plunge and also allows to determine the shift 
in the position of the binary during the inspiral. Our analytic approach is suitable for discussing 
the dependence of the recoil on various parameters characterizing the binary. 

In order to gain higher accuracy, the results derived in this paper can be further generalized 
by the inclusion of higher order corrections: the spin-spin, mass quadrupole - mass monopole, 
and second order relativistic contributions. 
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